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Gauge/gravity duality has provided unprecedented opportunities to study dy¬ 
namics in certain strongly coupled gauge theories. This review aims to highlight 
several applications to heavy ion collisions including far-from-equilibrium dynam¬ 
ics, hydrodynamics and jet energy loss at strong coupling. 

1. Introduction 

Gauge/gravity duality equates certain gauge theories with theories of quantum grav¬ 
ity in one higher dimension.^ This extra dimension has a natural interpretation as 
the renormalization group scale, and has led to the name “holographic duality." In 
the limit where the gauge theory is strongly coupled and has a large number of 
colors Nc the dual description reduces to classical supergravity. Hence, challenging 
strongly coupled quantum dynamics in the gauge theory can be accessed by solving 
classical partial differential equations. All physics — from microscopic interactions 
to macroscopic hydrodynamics — is encoded in the dual classical dynamics. Holog¬ 
raphy can thereby provide systematic and controlled access to strongly coupled 
dynamics via the solution to the classical partial differential equations. Currently 
there is no other theoretical tool capable of accessing strongly coupled real-time 
dynamics in in a controlled setting. 

Heavy ion collisions at RHIC and the LHC have demonstrated that the produced 
quark-gluon plasma is strongly coupled.^ Hence it is of interest to model heavy ion 
collisions using holography. However, while there exists many theories with dual 
gravitational descriptions — some with QCD-like features such as confinement and 
chiral symmetry breaking^’^ — the dual description of QCD is not known (if it 
exists). Nevertheless, given the immense challenge of studying strongly coupled 
dynamics in QCD it is invaluable to have a model where strongly coupled dynamics 
can be studied in a controlled systematic setting. This is especially true if the 
results are valid for a class of strongly coupled gauge theories and hence have some 
degree of universality. 
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In this review we shall focus on the simplest theory with a holographic dual: 3+1 
dimensional super-Yang-Mills theory with four supersymmetries {Af = 4 SYM). In 
the limit where the SYM is strongly coupled and has a large number of colors the 
dual description reduces to classical supergravity in asymptotically five dimensional 
Anti-de Sitter spacetime AdSs^ (see also®’®). While the ground state of SYM is very 
different from QCD — SYM is conformal and contains no particles whereas QCD is 
confining and has a rich spectrum of hadrons — at temperatures not too far above 
the deconfinement transition both theories consist of strongly coupled non-Abelian 
plasmas. It is in this setting where dynamics in SYM can potentially provide lessons 
for dynamics in QCD. 

In what follows we highlight insights obtained via holography in three settings. 
First, we discuss relativistic hydrodynamics including computations of the shear 
viscosity at strong coupling and the construction of causal viscous hydrodynamics. 
We then discuss far-from-equilibrium dynamics and the collisions of shock waves in 
SYM and their application to early time dynamics in heavy-ion collisions. Finally, 
we give a brief review of progress made in jet quenching and energy loss at strong 
coupling. 


2. First successes: viscosity and relativistic hydrodynamics 
2.1. Viscosity from black hole horizons 

Gauge/gravity duality provides a very natural way to describe thermal states: they 
are described by black holes, which are the simplest stationary states in AdS^ gravity 
with a notion of temperature. This temperature is simply the black hole’s Hawking 
temperature.^ The topology of the horizon is inherited from the gauge theory, 
where a thermal state usually is extended infinitely; this is one major difference 
from ordinary black holes in our universe, and in fact these black holes in AdS^ are 
more appropriately called black branes. Another major difference is the boundary 
of the AdS^ spacetime, which can be thought of as the place where the gauge theory 
lives. This boundary reflects back the thermal radiation, and the black branes are 
hence in stable thermal equilibrium, as opposed to black holes in flat spacetime, 
which evaporate away®'. 

It is also possible to construct gradient expansion solutions to Einstein’s equa¬ 
tions, where the local geometry is approximately that of boosted black brane with 
slowly varying temperature and boost velocity.®’® In this case the corresponding 
dual gauge theory state is that of a system in local thermal equilibrium with dy¬ 
namics governed by hydrodynamics. Simply put, in the limit of slowly varying 
fields Einstein’s equations reduce to hydrodynamics in one less dimension. The 

^Also, black holes in our universe tend to have very small temperatures of order of mK, as black 
holes of higher temperature are hard to form and would quickly evaporate. The thermal stability of 
black holes in AdS allows any temperature, and in fact the temperature of black holes representing 
a quark-gluon plasma would have the same very high temperature of order lO^^K. 
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connection between slowly varying solutions to Einstein’s equations and hydrody¬ 
namics is known as the fluid/gravity correspondence. By matching the gradient 
expansion obtained by solving Einstein’s equations to the constitutive relations of 
hydrodynamics, transport coefficients such as the viscosity can be computed at 
strong coupling. 

At first order in gradients the constitutive relation for the stress tensor in 
relativistic neutral hydrodynamics read^^ 


= eu^Uu -l-p[e]A^,, -I- where, (1) 

+ Ufj,Uy and ( 2 ) 

= -? 7 [e] - C[e] • u) + O(d^), with (3) 

- ?A^,A„^V“u^, (4) 


where e is the proper energy density, the fluid velocity, p[e] is the pressure 
(which is a function of the proper energy via an equation of state), is the shear 
tensor, 77 the shear viscosity and C, is the bulk viscosity, which vanishes in scale 
invariant theories such as SYM. Note that the fluid velocity and energy density are 
defined as the time-like eigenvector and associated eigenvalue of the stress tensor 
{Tfj.vu'' = —eu^) and that is transverse and traceless: = 0. 

Alternatively one can say that the fluid velocity is defined such that when boosting 
with velocity there is no momentum flow, i.e. Tq^ = 0, which is called the 
Landau frame'^. The conservation equation = 0 together with the equation 

of state and the viscosity (computed theoretically or measured) now form a closed 
system of equations, fully determined by a given initial energy density and fluid 
velocity. 

The equation of state, viscosity and higher order transport coefficients are not 
determined by hydrodynamics itself, and are a property of the microscopic theory 
under consideration. In weakly coupled QCD it is possible to use the Boltzmann 
equation to perturbatively compute the shear viscosity, which led to^^ 


y3 I 

^weak = K yrq r + ( ~~i ), 

r iog(i/5) r 


( 5 ) 


with T the temperature, g the gauge coupling and k a coefficient depending on 
the number of colors and flavors, which can be computed numerically. Non- 
perturbatively it is unknown how to compute the viscosity of QCD. However, there 
have been attempts to compute the viscosity using lattice techniques (see^^ for a 
review). It was therefore ground-breaking when a few years after the discovery of 
gauge/gravity duality Policastro, Son and Starinets computed the shear viscosity 
of a thermal plasma described by strongly coupled SYM theory: 

( 6 ) 


^When including a charge density with corresponding current and chemical potential it is also 
possible to define a frame where there is no charge flow, which is called the Eckart frame. In 
this review we do not include a charge density. 
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with s is the entropy density. The factor of 1 /47r is surprisingly universal and valid 
for any strong coupling and large Nf. quantum field theory with a gravitational dual. 
For instance adding charge density changes both the shear viscosity and the entropy 
density, but its ratio is still l/47r. 

Nevertheless, it is possible to have holographic models with different viscosity 
to entropy ratios, most notably by including higher derivative corrections in the 
Einstein equations. In this case corrections for finite can actually make the rj/s 
ratio drop below l/47r, thereby violating the conjectured universal lower bound^^ on 
the viscosity.^® Nevertheless, general considerations using the quantum mechanical 
uncertainty principle still hint at the existence of a lower bound, even though its 
value is likely smaller than l/47r. 

Holography can also yield finite coupling and finite Nc corrections to the viscos¬ 
ity. For SYM the results readd’’^^® 

1 / 1 c; c \i/2\ 

with A = g^Nc the’t Hooft coupling. As expected from (5) the viscosity increases 
as the coupling constant decreases, and the increase can be significant at reasonable 
values of the coupling constant (for A ~ 20 and Ac ~ 3 the corrections both are 
roughly 20%). Current hydrodynamic models of the quark-gluon plasma indeed 
favor an g/s in the range l/47r to 2.5/47r.^'^“^^ 

2.2. Second order relativistic hydrodynamics 




Figure 1. (left) The second order hydrodynamics derived using AdS/CFT (AdS) is compared 
with the older weakly coupled Muller-Israel-Stewart hydrodynamcis (IS), which uses Ttt = QijIsT 
and does not include the Ai or A 2 in eqn. 8. For these initial conditions the spatial and momentum 
anisotropy {ex and ep) do not depend much on the precise value of the second order transport 
coefficients (figure from^'^). (right) Similar hydrodynamics has been used to obtain recent viscosity 
over entropy ratios of roughly 0.12 and 0.2 at top RHIC and LHC collisions respectively (LHC 
estimate shown, figure from^^). 


Holography and in particular the fluid/gravity duality also has led to a more system¬ 
atic understanding of (relativistic) hydrodynamics itself.®’^^’^^ As alluded to above, 
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via the fluid/gravity correspondence it is possible to systematically solve Einstein’s 
equations. The gradient expansion solution encodes the constitutive relations of 
hydrodynamics to all orders in gradients and associated transport coefficients.®’^® 
Let us begin by simply stating the results of®’^® for the shear tensor up to the 
second order in gradients. 

Since SYM is conformal, to write at second order in gradients it is useful 
to employ derivative operators which transform covariantly under conformal 
transformations. For tensors Q^ ' which, under a conformal transformation —)■ 

transform like Q'i ' = T) is defined by^® 


Av = 


+ wAMi::. + + ■■■- iLQx::. 


A^ — u ditU^ 






With this definition derivatives transform covariantly: T>\Q^ - = -, 

and for instance Eq. (4) simplifies to + 'Di.u^. The full second order 

conformal shear tensor is now given by'’: 


^ ^ fiiy r g T-jy U nxfj-\- Ai[(J^^CT^ 


A 


z —cr 


'A 





w 


Ol^ 


^a/?] j 


( 8 ) 


where oj^i, = In strongly coupled SYM the second order transport 

coefficients are given by T,r = (2 — ln(2))/27rT, Ai = g/27:T = 2 A 2 /ln( 2 ) and 
A 3 = 0.® Note that all these coefficients get considerably more complicated when 
including a chemical potential; in particular A 3 becomes non-zero. 

Second order hydrodynamics is particularly useful for numerical simulations of 
viscous hydrodynamics. This is because high momentum modes in first order rela¬ 
tivistic hydrodynamics can propagate faster than the speed of light. Fundamentally 
this is not a problem, as hydrodynamics is an expansion in long wavelengths, i.e. 
small momenta, and these acausal modes are outside the regime of its applicabil¬ 
ity.^® Nevertheless, these modes make the equations numerically unstable^® and 
therefore in numerical simulations a second order term was added to stabilize the 
equations, giving Muller-Israel-Stewart hydrodynamics.®®’®® 

Muller-Israel-Stewart hydrodynamics can be seen as a special case of second 
order hydrodynamics where Ai = 0. One insight which holography helped solidify is 
that, in contrast to the Muller-Israel-Stewart theory, at a given order in the gradient 
expansion one must add every possible tensor structure to the constitutive relations 
which is consistent with the symmetries of the underlying quantum field theory. This 
had lead to the development of causal viscous relativistic hydrodynamics.^^’®^’®® In 

'^Most of these computations have been done for fluids living in a curved spacetime,^® which gives 
an extra term not considered in this review. 
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conformal theories such as SYM the equations of second order hydrodynamics give 
causal propagation for high momenta modes. 

While the second order transport coefficients make simulations much easier to 
perform, the precise value of these coefficients fortunately does not seem to have 
a big influence on the particle spectra in realistic simulations (see figure 1, left). 
The latter feature does however crucially depend on using conformal hydrodynamic 
equations; when using (second order) non-conformal equations results do seem to 
depend on the second order coefficients.^^ One of the biggest successes of these 
equations are estimates for the shear viscosity of the quark-gluon plasma,as 
shown in Fig. 1 (right). 

Another more recent example of a better understanding of hydrodynamics comes 
from,^® where in a specific boost-invariant setting all transport coefficients were 
computed numerically up to order 240 in the gradient expansion. With those coeffi¬ 
cients it could be shown that the hydrodynamic gradient expansion is not necessarily 
convergent, it was found to be an asymptotic series. Moreover, it was understood 
why hydrodynamics did not converge, since the divergent hydrodynamic expansion 
contains information about the non-hydrodynamic modes. Within AdS/CFT these 
modes are included in the so-called quasi-normal modes: vibrations of the black 
hole horizon. Strikingly, by a Borel resummation reference^^ was able to extract 
the precise value of the dominant quasi-normal mode from the hydrodynamic ex¬ 
pansion, thereby showing how hydrodynamics contains information about its own 
break-down. 

Much work has been devoted to study charged hydrodynamics, extending the 
previous formulas to a plasma with a conserved current, chemical potential and var¬ 
ious new transport coefficients, see for instance. More recently there has been 
interest in anomalous hydrodynamics, where quantum anomalies lead to modified 
hydrodynamic equations. A plasma with a non-trivial axial charge density will de¬ 
velop an electric current J in the direction of an external magnetic field B (chiral 
magnetic effect or CME):^"^’^® 

J= (9) 

where /is is the chemical potential associated with the axial charge. For averaged 
heavy ion events /is would be zero, but a full event-by-event anomalous hydrody¬ 
namic simulation suggests that this effect is measurable in heavy ion collisions.^® 
Similarly to the analysis above, also for the CME holography has proved extremely 
useful to show how this effect is in fact necessarily included already in first order 
charged hydrodynamics, and also to estimate the relevant transport coefficients at 
strong coupling. 

Lastly, inspiration from the symmetries present in AdS has led to non-trivial 
analytic solutions of the hydrodynamic equations. An interesting solution was 
provided by Gubser,^® where a solution to viscous conformal hydrodynamics with 
boost-invariance and a non-trivial expansion in the transverse plane was found. 
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This solution follows from a simple set-up in AdS, but then leads to this non-trivial 
analytic and hence very useful hydrodynamic solution. Recently, this approach has 
been extended to an exact solution of the relativistic Boltzmann equation. 

Interestingly, aside from the computation of transport coefficients, much of the 
development in hydrodynamics above need not rely on gauge/gravity duality. Nev¬ 
ertheless, inspiration from the duality has led to this more precise understanding of 
hydrodynamics, by using gauge/gravity as a theoretical arena to get new insights. 

3. Fast thermalization and the resulting stress tensor 
3.1. Simple models 

Perhaps an even more remarkable feature of holography is the straightforward pos¬ 
sibility to study far-from-equilibrium dynamics in strongly coupled gauge theories. 
How do non-equilibrium systems thermalize and relax to local equilibrium? Us¬ 
ing holography answering this question amounts to studying real time dynamics 
in the gravitational dual, which is not much different from conventional problems 
in general relativity, such as binary neutron stars, or black hole inspirals.As 
is well-known these latter problems can still be very challenging from a numerical 
point of view, but compared to non-perturbative quantum field theory it is defi¬ 
nitely an enormous simplification. For a review of the use of numerical relativity in 
holography see for example.^® 

The first studies of holographic thermalization computed how perturbed black 
holes relax to homogeneous black holes, computing the so-called quasi-normal 
modes.These are determined by expanding the perturbation of the metric in 
a sum of modes proportional to with complex frequency uji. The Einstein 

equations and boundary conditions at the horizon and boundary then dictate that 
there is a discrete spectrum, whereby the mode always decays (i.e. $ 5 ( 0 ;^) < 0). The 
mode with the slowest decay will then determine the late-time decay back to the 
thermal black hole, and hence can be said to determine the thermalization time in 
this setting. Dimensional analysis in this scale invariant theory dictates that Wi ~ T 
, but the coefficient has to be computed numerically, giving 9(wo) ~ 8.6T.^^ 

The estimate from quasi-normal modes could imply a thermalization time 
ttherm < 1/T, which for a temperature of 1000 MeV would be about 0.2 fm/c. 
Compared to perturbative methods this is very fast, which should be expected for a 
system at strong coupling. On the other hand, even though the coupling is assumed 
to be infinite, there is still a finite answer for the thermalization time, leading to 
the intuition that strong interactions still have some microscopic physics present as 
the scale 1/T. 

Nevertheless, the the quasi-normal mode analysis is only valid near equilibrium: 
the state is assumed to be a small perturbation of the final thermal state. In heavy 
ion collisions the initial state is very far-from-equilibrium, so solving the full non¬ 
linear Einstein equations is of great importance. One of the first studies solving this 
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numerically is presented in,^®’"^® where quark-gluon plasma formation and thermal- 
ization was studied in both homogeneous and boost invariant set-ups (see fig. 2). 
The gauge theory dynamics was, by symmetry, effectively 1-1-0 dimensional and the 
dual gravitational dynamics was 1-1-1 dimensional. In these papers non-equilibrium 
states were created by starting in the ground state and turning on background fields. 
It was found that even in the case where the initial state is the vacuum, a thermal 
state can be achieved in times of order 1/T with T the final state temperature. 

An intriguing finding in"^® was that when the system first starts to behave hydro- 
dynamically, viscous effects were order one: the first order viscous correction to the 
constitutive relations (1) is as big as the leading term. Nevertheless, even with such 
large viscous corrections the hydrodynamic gradient expansion was seen to be well 
behaved. Moreover, the applicability of hydrodynamics was seen not to be governed 
by convergence of the hydrodynamic gradient expansion. Instead, the applicabil¬ 
ity of hydrodynamics was seen to be governed by the decay of non-hydrodynamic 
degrees of freedom. 

The approach in^® was extended in.®'^ Here the main feature was the possibility 
to chose several initial states and to evolve them without turning on any background 
fields. For all initial states they found a thermalization time less than 1 /T, with T 
the temperature at the moment of thermalization. 

Of particular utility, this paper also introduced a dimensionless quantitiy, 
F(w) = ^drW, where w = T^ffT can be thought of as dimensionless time and 
Te/y is the temperature associated with the energy density at that time. The di¬ 
mensionless function F{w) is completely determined within hydrodynamics and for 
instance equals 2/3 for ideal hydrodynamics. As F{w) is unique within hydrody¬ 
namics it is clear that differences between F(w) found for the different evolutions 
in AdS/CFT must be due to genuine non-hydrodynamic effects, as is convincingly 
shown on the left of Fig. 3. 

The right of figure 3 shows a similar study in a homogeneous setting, with more 
than 1000 initial states, also leading to thermalization within a time of 1/T.®^’®^ 
The main new result in this paper was to compare the full non-linear thermaliza¬ 
tion process with the linearized quasi-normal mode approach outlined in the second 
paragraph of this section. Quite surprisingly it was found that for general initial 
states the evolution was always well described, within 20%, by the linearized ap¬ 
proximation. This then also confirmed that the quasi-normal mode estimates of 
strongly coupled thermalization times were indeed accurate also at the non-linear 
level. 

Lastly, these studies can be corrected for finite coupling effects, which has re¬ 
cently been done in.^®’®^ Not surprisingly, finite coupling corrections slow down the 
thermalization, as in these papers becomes apparent through a smaller imaginary 
part of the lowest quasi-normal mode. 
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Z 


Figure 2. A simple and often used model of a heavy-ion collision was proposed by Bjorken,^^ 
where he assumed that a heavy-ion collision is approximately boost-invariant, at least near 2 = 0. 
This means that all physics only depends on proper time r = y/t'^ — 2 ^, and hence that all physics 
experienced by observers, such as a temperature field illustrated here, is independent of its frame. 
The two red lines would illustrate two such frames, which indeed has equal temperatures at equal 
r (figure from^®). 



Figure 3. The left figure shows boost invariant thermalization, where the evolution of 29 random 
initial states of F{w) = as a function of it; = T^ffr is compared with predictions from first-, 

second- and third-order hydrodynamics (figure from^^). It can be seen that they all thermalize in 
a time where r < O.^/T^ff, where T^ff is the temperature associated with the energy density at 
that time. The right figure did a similar analysis in a homogeneous set-up, using over 1000 initial 
states. Also there thermalization always happens within t < 1/T. Here it is also shown that the 
difference between the linearized and fully non-linear Einstein equations is small, almost always 
less than 20% (figure from^^). 

3.2. Colliding shock waves 

All models presented above contained enough spacetime symmetry to make the 
gauge theory dynamics effectively 1 + 0 dimensional. Here we shall review results 
where the gauge theory dynamics are allowed to be 1 + 1 and 1 + 2 dimensional. In 
holography the energy-momentum in the gauge theory directly corresponds to the 
metric in the gravitational theory, so to simulate highly energetic colliding nuclei 
it is therefore natural to consider collisions of gravitational shock waves, moving at 
the speed of light.The shock waves can be viewed as analogous to the shock 
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waves constructed in flat space by Dray and ’t Hooft,®^ which consist of taking a 
Schwarzschild black hole and boosting this to the speed of light, keeping its total 
energy fixed. The shock waves are dual to lumps of energy in the gauge theory, 
which propagate at the speed of light. When the shock waves collide they can 
produce a black hole which will then relax to equilibrium. Likewise, in the dual 
gauge theory, the collision will result in the formation of a non-Abelian plasma 
which will subsequently thermalize. 

One interesting result from the aforementioned studies of shock wave collisions 
was the entropy production, estimated by the area of the trapped surface of the 
formed black hole at the time of collision. This trapped surface gives a lower bound 
on the area of the black hole horizon, which is directly related to the produced 
entropy in the gauge theory. Ref.®® estimated that: 

S ^ Strapped ~ 

where S is the entropy produced in a central gold nucleus-nucleus collision with a 
nucleon-nucleon collision of energy and Strapped is the entropy of the trapped 

surface. The entropy increase during the evolution of the quark-gluon plasma is 
not too big due too to the relatively small viscosity, but can in fact be sizable as 
highly anisotropic plasmas increase the entropy production. Nevertheless Eq. (10) 
gives a lower bound on the final entropy produced, which is directly related to the 
measurable total charged particle multiplicity Ntot,ch = -S'/?.5.®^ In®® the estimate 
(10) was made somewhat stronger by doing calculations numerically till later times, 
but two important features remained. Firstly the numerical value is about right 
for RHIC collisions, albeit on the high side, which is impressive as Eq. (10) is a 
direct calculation, without having any input from other experiments. Secondly, 
the 2/3 power law as a function of y/sj^jy is a robust outcome of a scale invariant 
theory. Nevertheless, with the LHC results now known,®® the experimental data 
favors a power of 1/2, and Eq. (10) therefore overestimated the multiplicity at LHC 
collisions. This can be taken as an indication that QCD is not scale invariant at 
energies probed between RHIC and LHC. 

The previous calculations could not say much about the dynamics after the colli¬ 
sion, which requires solving the full non-linear Einstein equations. It was therefore a 
major breakthrough when a real collision could be simulated within gravity.®® This 
study collides planar gravitational shock waves, thereby simplifying the problem 
by neglecting transverse dynamics. The initial conditions for the shock wave are 
determined by the stress tensor 


T±±{z±) = Kh{z±), (11) 

where z± = t ± z are null coordinates along the collision axis and k = fV//27r^ is 
a measure of the degrees of freedom of the gauge theory. All other components of 
the stress tensor are zero, and within pure gravity this fully determines the metric 
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before the collision. For this review we restrict h to be Gaussian: 


h{z) = 




exp 


2 " 


( 12 ) 


where w is the width of the Gaussian and Kfj,^ = K h{z)dz is the total energy 
per transverse area. 

In what follows we make the following choices for parameters. We choose W = 
1.8 such that the ratio e/T^ « 12 matches lattice QGD simulations.^® The choice of 
W < 3 compensates for the fact that SYM has approximately three times as many 
degrees of freedom than QGD. We then fix /i by setting the 




y/SNN Ypb 


(13) 


where IVpb = 206 and R = 6 fm. The right hand side of this equation is simply 
the energy density per unit area near the center of a lead nucleus. We choose the 
width of the Gaussian to be w « | Rtun/, with niN the nucleon mass. Of 
course real nuclei are not Gaussians, (nor is QGD = SYM) so all these numbers 
have to be taken with a grain of salt. Nevertheless, these choices gives a qualitative 
understanding of the scales in the problem. 

After setting up these initial conditions it is possible to numerically solve the 
Einstein equations (see^®’®® for an introduction), and thereafter extract the expec¬ 
tation value of the gauge theory stress tensor. Figure 4 shows the energy density 
obtained in this way for collisions having = 19.3 GeV and ^smn = 2.76 TeV. 

For a proper comparison with®®’®^ one should notice that the dimensionless product 

1/3 

= (^ ^2/1^ ttR ^ ) 3 R^n/V^nn ~ 23(-\/sArAr/GeV)“^^®, (14) 


so that 19.3 GeV and 2.76 TeV correspond to a of 3.2 and 0.11 respectively, 
whereby we stress once more that these are ballpark figures, which we think are 
helpful for a useful comparison of gauge/gravity results to actual heavy ion collisions. 

Several curious features are visible in the resulting energy density.®^ Firstly, for 
the low energy collisions the peak energy density is about 60% higher than the sum 
of the incoming energies. This is interpreted as a ‘pile-up’ of energy, as the shocks 
thermalize during the collision. For high energy shocks there is not enough time for 
this effect, and the shocks therefore only lose there energy gradually to the plasma, 
as well as forming a transient region of negative energy density, where it is not even 
possible to define a local rest frame.®® Guriously, for the high energy shocks there 
seems to be a reasonably accurate analytic formulation — at least at small rapidities 
— by using a complexified boost symmetry, even valid in the far-from-equilibrium 
regime.®® At later times the low and high energy collisions behave similarly though; 
they both thermalize fast (shown in figure 5) and deposit all their energy into the 
plasma. 

Interestingly, the formed plasma is quite insensitive to the profile of the initial 
shocks, provided characteristic length of the profile is not too big.^® This makes 
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the local energy density of the resulting plasma robust at least for thin and hence 
high energy shocks. The rapidity profile of these collisions may even be universal 
at strong coupling, and is therefore of considerable interest (figure 6, right). The 
left part of that figure displays the similar profile for low energy collisions, which 
looks similar in shape, but with a width that becomes smaller as the collision energy 
decreases. 



Figure 4. Energy density e for collisions of shock waves at low energy (.y/s^vv = 19-3 GeV, left) 
and high energy (y/sjvlv = 2.76 TeV, right) as a function of time and longitudinal coordinate z. 
The grey planes lie at the origin of the vertical axes. Note that these simulations are planar in the 
transverse plane and therefore neglect transverse dynamics. For the short time scales presented 
at the right plot this is probably a good approximation, but the simulation on the left will change 
when taking transverse expansion into account as well (figures adapted from®^). 



Figure 5. 3AP)°‘^/£lioc for low (left) and high (right) energy shocks, with AP the difference 

between the longitudinal pressure from holography, and the one computed using eqn. 1. The 
white areas indicate regions outside the light cone or where hydrodynamics deviates by more than 
100%. The black lines are again stream lines, which for the high energy shocks stop at the region 
where no local rest frame exists. 
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Figure 6. Energy density in the local rest frame around mid-rapidity as a function of spacetime 
rapidity y and proper time r for thick (left) and thin (right) shocks (the same shocks as figure 4). 
In the right case we have excluded from the plot the region in which the local rest frame is not 
defined because 2|s| > |e-|- Pl\ (figures adapted from®^). 

3.3. Including transverse dynamics 

Lastly, we would like to mention that it is also possible to include transverse dynam¬ 
ics. Transverse expansion is essential during the late time of the QGPs explosion, 
but also early time initial velocities in the transverse plane influence the transverse 
momenta spectra. Nevertheless, a full shock wave collision with with transverse 
expansion requires a 3+1 dimensional code in AdS spacetime, as there is always the 
extra AdS direction. As this can be technically complicated it is therefore natural 
to first try to include transverse dependence while assuming boost invariance in the 
longitudinal direction. 

Assuming boost invariance is a bit subtle, as Fig. 6 convincingly shows a signif¬ 
icant dependence on rapidity. Also, using proper time as a variable usually neces¬ 
sitates to start the gravitational simulation a finite time after the collision (see^^ 
for a counter example though), which has to be accompanied by an assumption on 
the dynamics before this time (early time expansion in the case to be presented). 
Nevertheless, starting the simulation at such a finite time with several simple pro¬ 
files gave two lessons.Firstly the plasma again thermalized fast, at a similar 
time scale as the examples above. Secondly, after the thermalization the transverse 
velocity was given to a good approximation by the following approximation:^^’^^ 


Vi = —0.33r ditje, 


(15) 


with T the proper time, and e the local energy density in the transverse plane. This 
formula is in spirit similar to the universal pre-flow found in,^® but note that only 
the flux is universal in the sense of that paper, whereas here we are interested in 
the local fluid velocity, which depends on the pressures as well. 

Using these simulations it was then possible to couple this to relatively standard 
hydrodynamic and cascade codes, and thereafter obtain the full transverse momenta 
spectra at mid-rapidity for a central collision, which fits data surprisingly well (figure 
7). Recently this has also been extended to several other systems, most notably 
also computing the measurable HBT radii. 
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Figure 7. From the local energy density we plot the temperature (assuming equilibrium, therefore 
labeled Tpseudo) and the radial velocity Vp = Upjur for a representative simulation. The plot 
illustrates four physical tools used: i) early time expansion, ii) numerical AdS evolution, iii) viscous 
hydrodynamics until T = 0.17 GeV, iv) kinetic theory after conversion into particles (indicated by 
arrows). The right plot shows the resulting particle spectra for pions and kaons, and also includes 
curves for a zero flow model (ZF) with Vi = 0 until t = 1.0 fm/c, and a free streaming model (FS), 
with zero longitudinal pressure (figure from,^^ experimental data from ALICE^®). 


4. Jet energy loss 

4.1. Are there jets at strong coupling? 

During the early stages of heavy-ion collisions — in which the energy density is the 
highest — energetic partons can be created via hard processes. The resulting exci¬ 
tations — jets — consist of collimated sprays of energy and can traverse the fireball 
while depositing energy and momentum into the medium. Analysis of the particle 
correlations in the produced jets can provide useful information about the dynamics 
of the plasma including the rates of energy loss and momentum broadening. 

Why study jets and energy loss at strong coupling? One piece of evidence that 
the quark-gluon plasma produced in heavy-ion collision is strongly coupled comes 
from jet quenching itself.®^^®^ Jets traversing the fireball appear to lose their energy 
very quickly, indicating strong interactions with the medium. 

A natural process by which jets can be created is the decay of off-shell particles 
such as photons. At weak coupling the decay leads to well-collimated sprays of 
energy with the polarization of the photon imprinting itself in the final state in a 
so called “antenna pattern" (see for example®"^ and references therein). However, at 
strong coupling the situation is very different: the decay of off-shell photons leads 
to states which are spherically symmetric (in the rest frame of the off-shell photon): 
all correlation functions at spatial infinity are spherically symmetric and there are 
no well collimated sprays of energy.®® 

What does this mean for studying jets at strong coupling? Does the observed 
isotropy mean there are no jets at strong coupling? Where does the isotropy come 
from? In other words, at what stage in the evolution does the isotropy develop? 
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Does it come from the underlying photon decay or from the subsequent late time 
evolution as the energy propagates outward to infinity? 

In order to understand the above issues it was useful to study the propagation of 
radiation in strongly coupled holographic gauge theories. A simple setting to study 
this problem is that of radiation produced by an accelerated quark in strongly 
coupled conformal field theories. If one takes a heavy quark and accelerates it along 
some given trajectory, what does the pattern of emitted radiation look like? This 
problem is analogous to textbook studies of radiation produced by an accelerated 
charge in classical electrodynamics. However, instead of radiation carried by a 
classical U{1) gauge field, the radiation is that of a strongly coupled non-Abelian 
quantum field theory. Is the emitted radiation isotropic at strong coupling? If not, 
does it isotropize as it propagates to infinity? 

This problem was studied in.®®’®"^ There it was found that the radiation produced 
by the accelerated quark is not isotropic. Nor does the radiation isotropize as it 
propagates out to infinity. Remarkably, the angular distribution of power produced 
by an accelerated quark is the same at weak and strong coupling! This suggests 
that the isotropy observed in the decay of off-shell photons at strong coupling comes 
from the production mechanism itself — i.e. the decay — and not the subsequent 
evolution. 

The lesson from this is that one should not treat the jet production mechanism 
with strong coupling tools. Indeed, in asymptotically free QCD it is natural to 
expect the production mechanism during the early stages of heavy ion collisions to 
be governed by weakly coupled physics. Ideally, one should treat the initial hard 
physics with perturbative QCD and then study soft processes such as energy loss 
with strongly coupled tools. Simply put, the utility of modeling jets with strongly 
coupled tools comes from studying their evolution and energy loss after they are 
produced. 

4.2. Holographic models of jets 

While at strong coupling it is impossible to create well collimated sprays of energy 
via off-shell photon decays, one can construct states at some initial time t = U 
which share features with jets created via weakly coupled processes. Namely, one 
can construct states with excitations which are well-collimated, have energy E ^ T 
with T the temperature of the plasma, and which can propagate very far through 
the plasma before thermalizing. Ideally, one would compute the initial state using 
perturbative QCD. However, we shall take a different approach and focus on univer¬ 
sal features which are insensitive to the precise initial conditions used to construct 
the state at time U. These include the rates of energy loss and momentum broad¬ 
ening and the penetration depth defined as the total distance a jet can propagate 
through plasma before thermalizing. 

Let us focus for simplicity on qq jets in strongly coupled conformal field the¬ 
ories (for other holographic models of jets see for example®®^®^). According to 
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holographic duality a qq pair moving through plasma is dual to a classical string 
moving through a black hole geometry.^ The geometry dual to an infinite static 
plasma at temperature T is the AdS-Schwarzschild geometry, where the metric may 
be written 


ds^ 




—f{u)dt^ + dx^ 


dv? 

fiu) 


(16) 


where L is the AdS radius and f{u) = 1 — u^lu\. The horizon of the geometry is at 
radial coordinate u = Uh with Uh = ^ /ttT. The boundary of the geometry, which is 
where the dual field theory lives, is at radial coordinate u = 0. 

As it evolves in time the string sweeps out and 1 + 1 dimensional worldsheet. 
Events on the worldsheet can be labeled by coordinates (r, a) which are related to 
events in spacetime via the embedding functions = {t{T,a),x{T,a),u{T,a)}. 
The dynamics of the embedding functions are governed by the Nambu-Goto action 


S = -To 


drday^—'f, 


(17) 


where the string tension is Tg = with A the’t Hooft coupling and 7 = det 'jab 
with 


7a6 = daX ■ dbX, (18) 

the worldsheet metric. Here and below the indices (a, b) run over the worldsheet 
coordinates r and a. The action (17) is simply the area of the worldsheet. 

Varying the action we obtain the string equations of motion 

drUi + d^u-, (19) 

where the string worldsheet currents H^ are 

= s{dlxM) = -ToV^r'^dMNdaX^, (20) 

where gMN is the AdS-Schwarzschild metric and as usual 7 “^ is the inverse of ^ab- 
Note that the energy of the string is 

E = -JdaUJ^. (21) 

The string equations of motion (19) are simply the equations of energy conservation 
on the worldsheet with Hq the energy flux down the string. 

The endpoints of the string are charged under a bulk C/(l) gauge field Am- The 
boundary of the AdS geometry behaves as an ideal conductor.®^ Hence the gauge 
field induces a current on the boundary, which up to an overall sign has the inter¬ 
pretation as the baryon current of a dressed quark.This is shown schematically 
in Fig. 8 . Therefore, one can regard the spatial location of the string endpoint as 
roughly coinciding with the spatial location of the quarks. The boundary conditions 
at the string endpoints are set by the mass of the quark.For infinitely massive 
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X 


Figure 8. A cartoon showing a string at some fixed time in the AdS-Schwarzschild geometry. 
The endpoint of the string is charged under a f/(l) gauge field Am- The boundary of the AdS- 
Schwarzschild geometry — located at radial coordinate u = 0 — behaves as an electromagnetic 
conductor. Up to a minus sign the current induced on the boundary has the physical interpre¬ 
tation as the baryon current of the quark. Likewise, via Einstein’s equations the string perturbs 
the AdS-Schwarzschild geometry. The near-boundary behavior of the metric perturbation encodes 
the stress tensor of the jet. 


quarks the endpoint is fixed at the boundary. Correspondingly, the current is 
localized with delta function support. In contrast, for massless quarks the endpoints 
are free to fall in the AdS-Schwarzschild geometry. Likewise, the current spreads 
out as the endpoints fall. For intermediate mass quarks the endpoint is fixed to be 
at some radial coordinate Um- 

Via Einstein’s equations the presence of the string also perturbs the AdS- 
Schwarzschild geometry. Just as the near boundary behavior of the gauge field 
Am encodes the quark baryon current J^, the near-boundary behavior of the met¬ 
ric perturbation Sqmn encodes the stress tensor of the system. This is also 
shown schematically in Fig. 8. 

With the above preliminaries layed out, the basic strategy then will be to con¬ 
struct initial string data at some time ti and evolve the string profile forward in 
time according to the string equations of motion (19). To echo the previous point, 
we seek strings whose boundary boundary energy density is that of a well colli¬ 
mated spray of energy which propagates very far before thermalizing. We will then 
study jet quenching and other observables which are insensitive to the details of 
the precise form of the initial data. In what follows we will highlight results for 
infinitely massive quark jets and massless quark jets. 











21st January 2015 1:22 


World Scientific Review Volume - 9.75in x 6.Sin 


chapterl.4 


18 


Paul M. Chester and Wilke van der Schee 


4.3. Heavy quark energy loss 


One of the first studies®^^®® of energy loss via AdS/CFT came from the study of 
heavy quarks propagating at constant velocity v through a static box of plasma 
at temperature T. As the quark moves through the plasma it transfers energy and 
momentum to the plasma via frictional drag. Hence an external force is required 
to maintain constant velocity. This can be supplied by an electromagnetic field 
coupled to the quark’s baryon current J^. 

Under the assumption that the quark has been dragged at constant velocity for 
an arbitrary duration of time, it is natural to seek a steady-state solution to the 
string equations of motion. Choosing worldsheet coordinates t = t and a = u the 
steady-state ansatz take the form 

x{t,u) = vt + vxq{u), (22) 


with the endpoint terminating at the boundary u = 0. Solving the string equations 
of motion (19) for xo(u) yields 


xo{u) 


1 

‘2uh 


-1 

tan — 
Uh 


+ 2 log 


Uh-U 
Uh + u 


(23) 


We plot xo(u) in Fig. 9. 

The electromagnetic field dragging the quark (and hence the string endpoint) 
must be supplying energy at a rate equal to the rate energy is dissipated by drag. 
In the dual gravitational description the energy lost via drag is encoded by a flux 
of energy down the string towards the event horizon, ^ = Hg. Plugging the the 
solution (23) into Eq. (20), the drag reads®'*^®® 


dE 

dt 


Hence, as the ’t Hooft coupling A —>■ 


'VXT'^ 


2 y/1 — V 

oo the quark loses energy faster and faster. 


(24) 



ttTx 



Figure 9. Left: the string profile (23) for a heavy quark. Right: the energy density produced 
by a heavy quark moving at three quarters the speed of light in the fc—direction. A Mach cone is 
clearly present. 


A natural question to ask is where does the energy lost by the quark go? To 
answer this question it is useful to study the stress tensor . This requires solving 
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Einstein’s equations for the metric perturbation SgMN due to the presence of the 
string and extracting the stress tensor from its near-boundary asymptotics. This 
calculation was carried out In the right panel of Fig. 9 we also plot the 

normalized energy density E = for a heavy quark moving in the x—direction 

at three quarters the speed of light. The speed of sound in a relativistic conformal 
field theory is Cg = l/-\/3 « 0.57. Hence the quark is moving supersonically. As is 
clear from the figure, the quark creates a Mach cone. Evidently, the energy lost via 
drag is carried away by sound waves. 

There have since been numerous studies of heavy quark energy loss at strong 
coupling in more general settings. Some examples include heavy quark energy loss 
in non-conformal theories,anisotropic plasmas^'^^ and in far-from-equilibrium 
states such as the shock wave collisions discussed above.In the latter study it was 
demonstrated that even in far-from-equilibrium states, such as those produced by 
colliding shock waves, the equilibrium energy loss formula (24) provides a reasonable 
estimate of the heavy quark energy loss rate. Additionally, heavy quark energy loss 
due to deceleration {i.e. non steady-state) was studied in,^'^"' and the rates of heavy 
quark momentum broadening were studied 

5. Light quark energy loss 

In contrast to strings dual to heavy quarks, where the endpoints are fixed to the 
AdS boundary, the endpoints of strings dual to massless quarks are allowed to fall 
unimpeded towards the event horizon of the AdS-Schwarzschild geometry.For 
massless quarks the string equations of motion (19) must be solved subject to the 
open string boundary conditions 

= 0 at the endpoints. (25) 

The fact that string endpoints fall toward the horizon implies there does not exist 
steady-state solutions to the string equations of motion. Hence studying light quark 
energy loss is considerably more complicated than studying heavy quark energy loss. 
Nevertheless it is a straightforward matter to numerically solve the string equations 
of motion (19) with the open string boundary conditions (25). For details on how 
to do this see for example.®® 

In Fig. 10 we plot a numerically generated string profile at several values of fixed 
coordinate time t. At time t = 0 the initial string profile is simply a point located at 
the AdS boundary u = 0. As time progresses the string expands from a point to an 
extended object and the endpoints fall towards the horizon. The endpoints travel 
a total distance Xstop ~ IO/ttT before falling into the horizon. As the endpoint fall 
into the horizon, the boundary baryon current (depicted schematically in Fig. 8) 
becomes delocalized and starts to evolve hydrodynamical according to the diffusion 
equation.®® This corresponds to the thermalization of the light quark jet. 

The total distance the string endpoints and hence the jet can travel depends on 
initial conditions. A natural question is for a fixed energy E, what is the maximum 
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ttTx 

Figure 10. A numerically generated falling string profile dual to a light quark jet. The string is 
shown at several values of fixed coordinate time t. The string is created at the point x = u = 0 
and subsequently expands into an extended object. The red lines denote the string endpoint 
trajectories. The endpoints travel a distance Xatop ~ IO/ttT before falling into the horizon at 
U = Uh. 


distance the string endpoints can travel before falling into the horizon? In the 
language of the dual field theory, for fixed jet energy E, what is the maximum 
distance a light quark jet can travel before thermalizing? To answer this question 
one can numerically generate an ensemble of solutions to the string equations of 
motion (19) with different initial conditions. One can then compute the energy (21) 
and compare the stopping distance aistop to E- This computation was carried out 
in®^ where it was found that the maximum distance a jet with energy E can travel 
before thermalizing is given by 

- f (tyl) ■ 


where C « 0.526. The Xgtop ~ E'^/^ scaling was also found in.^°^ The constant C 
was recently computed analytically and found to be^°® 


. 2i/3r(|) 

v^r(l) • 


(27) 


We now turn to the light quark energy loss rate. Following,we focus on strings 
whose Xstop is asymptotically large compared to 1/T. Such string have worldsheets 
which are approximately null.®^’^°® Why? When Xstop —>■ oo the scaling (26) 
requires E —>■ oo. Since strings have finite tension the E ^ oo limit is generically 
realized by strings that expand at nearly the speed of light, meaning that the string 
profile must be approximately that of an expanding filament of null dust. Indeed, 
null strings have profiles satisfying 7 (Xnuii) = 0 and from (20) have divergent 
energy density. 

Since null strings satisfy 7 (?'fnuii) = 0 they minimize the Nambu-Goto action 
(17) and are exact albeit singular solutions to the string equations of motion (19). 
Following, to obtain finite energy solutions to the equations of motion one can 
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expand the string embedding functions about a null string solution 

= XZ, + + e^SX^^ +..., (28) 

where e is a bookkeeping parameter (which can be set equal to 1 at the end of 
calculations). In what follows it is useful to choose worldsheet coordinates t = t 
and a such that Xnuii • -^nuii = ^ and = {0, 0}- We denote the location of 

the string endpoints by ct = cr*. The string equations of motion can then be solved 
perturbatively in powers of e. 

Focusing on strings propagating in the a;—direction, the null string embedding 
functions can be written 


^null — {^7 2^geo(t, O'), 0, 0, Ugeo(t, O')}, 


(29) 


where for each a Xgeo and rtgeo satisfy the null geodesic equations, which read 


_ J 

dt ~ C’ 

dUgeo ^ /\/C^ - / 

dt i 


(30) 

(31) 


where ^ = ^(cr). The parameter ^ determines the initial inclination of the geodesics 
in the x — u plane and, more fundamentally, specifies the conserved spatial momen¬ 
tum associated with the geodesics, f{u)~^dxgeo/dt = The geodesic equations 
have the solution 

^geo = -^2F^ (i, i; I; + Cia) (32) 

^eeo ' ^ 


where 2 F 1 is the Gauss hypergeometric function, ^ = 1/(1 — /'^) and C'(cr) is an 
arbitrary function which can be chosen such that Xgeo{t = 0, u) = 0. 

Henceforth let us focus on strings created at the point a: = t6 = t = 0at the 
AdS boundary. At asymptotically early times t uu, when the string is close to 
the AdS boundary, / = 1 and the geodesics are given by 


Xgeo = t cos cr, Ugeo = t sin O'. (33) 

Hence, the worldsheet coordinate cr is simply the initial angle of the geodesic in the 
X — u plane. Likewise, 


^(cr) = sec cr. 


(34) 


In Fig. 11 we plot a null string generated by a congruence of geodesics with 
cr* = 0.025. The string profile, denoted by the blue curves, is shown at several 
values of coordinate time t. The string starts off at a point on the boundary and 
expands at the speed of light while falling towards the horizon. The rust colored 
curves represent the null geodesics that make up the string and the red curve is the 
endpoint trajectory. The total distance the string endpoint can travel is entirely 
controlled by the angle cr*. As we elaborate on further below, in the limit cr* —)■ 0 
the stopping distance diverges. 
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Figure 11. A null string (blue) shown at several different coordinate times t. The string starts 
off at a point on the boundary and expands at the speed of light while falling towards the horizon. 
The rust colored curves represent the null geodesics that make up the string and the red curve is 
the endpoint trajectory. 


To study energy loss rate we must compute the first order corrections to the 
string embedding functions. At leading order in the bookkeeping parameter e the 
worldsheet energy density IIJ and flux IIq read 


K = - 


To^daU^ 


geo 




2efdtSx(^i) 


+ 0{y/e), IIq — 0{y/e). 


(35) 


Hence at leading order in e the string equations of motion (19) simply read cltHg = 0 
so nj = n5((T). In other words, energy is simply transported on the congruence of 
geodesics which make up the null string! 

We define the energy 

ranix) 

E{x) = - / n0(a), (36) 

J CT* 


where a nix) is the a corresponding the the geodesic which impacts the horizon at 
X. In particular, <jh{x = Xstop) = cr*. The energy E{x) is simply the total string 
energy that makes it past the point x. The energy loss rate per unit length is 

^ = -n5(crff(a:))cr^(a;). (37) 

We discuss the precise boundary interpretation of dE/dx below. To compute dE/dx 
we must compute Hq and anix). 

Because Hq is time independent it is only necessary to compute it at an initial 
time. Near the AdS boundary, where / = 1 and the geodesics are given by (33), 
the string equation of motion dtU^ = 0 leads to the equation of motion for Sx(i) 

2 2 

df fe(i) + jdtSx(i') = 0, 


( 38 ) 
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which has the solution 

dx{t, a) = (j){a) + (39) 

for arbitrary functions (j){(j) and 'ip{cr). The open string boundary conditions (25) 
require = 0. 

Substituting the string solution (39) and the geodesic solution (33) and (34) into 
(35) we obtain 


IIq = —To CSC (Ti 


/ CSC 2a sin a 


e-ipia) 


+ 0(v^)- 


(40) 


With nj computed and the congruence of null geodesics specified we can now 
compute dE/dx. But what is V’(f^)! The energy loss rate depends on an arbitrary 
function. At first sight this seems like a disaster. Not only can we not describe the 
production mechanism of jets via holography but their energy loss rate, at least for 
massless jets, seems to be very sensitive to initial conditions. However, this grim 
diagnosis turns out to be incorrect: the ambiguities in the light quark energy loss 
rate are largely transient effects which die out over time scales of order 1/T. These 
effects are negligible in the high energy limit when Xgtop oo. 

Let us consider Xgtop oo limit with x/xsto-p fixed and first compute anix). 
anix) can be computed from the geodesics solution (32) by setting Ugeo = Uh 
and solving the resulting transcendental equation for a. With ^ given by (34) the 
resulting equation can be solved analytically in the limit ct —)■ 0 (i. e. for geodesics 
which travel very far before falling into the horizon). The solution reads a nix) = 

■ Likewise, setting anix = a;stop) = cr*, we see that the stopping distance is 
related to cr* via 


_ «.r(i)2 

As previously advertised, Xgtop oo when cr* —>■ 0. We therefore may write 



(41) 


(42) 


Therefore, when we take Xgtop oo with a:/a;stop fixed, we are forced to take 
an —> 0. Simply put, when Xstop —t oo, energy is transported along geodesics 
which originate asymptotically close to the string endpoint. Only these geodesics 
propagate to distances x = O(xstop)- This behavior can clearly be seen in Fig. 10. 

Therefore, to get the energy loss rate we may expand Eq. (37) about an ~ cr* —> 
0. The function ip{a) may be approximated by ip^a) ss ip'{as,){a — (t»). Hence, near 
the string endpoint 


H 


0 


To 


a‘^yj2€'4)'{a^){a — ct *) 


(43) 
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The energy loss will depend on the unknown parameter ip'{as,) which can be repack¬ 
aged in terms of the initial energy 


E{x = 0) = 


TrTn 


+ 0(\/^/oPP). 


Substituting (43) and (41) into (37) and dividing by (44) we then secure^°® 

1 dE 


E{x = 0) dx 




(44) 


(45) 


Eq. (45) together with the stopping distance (26) provides a complete description 
of energy loss in terms of the initial energy, ’t Hooft coupling A and temperature T. 

We now return to the boundary interpretation of dE/dx. It turns out that it 
is possible to solve the gravitational backreaction problem analytically in the long 
wavelength limit.In doing so the hydrodynamic limit of the stress tensor 
can be extracted from the gravitational field perturbation in closed form. How is 
the string energy loss rate (45) encoded in The answer is 

^M'^hydro = - xt), (46) 

with = j^9{xstop — x). Simply put, the energy loss rate (45) describes the 

rate in which energy is transferred from the jet to hydrodynamic modes. 

The light quark energy loss rate (45) stands in stark contrast to the heavy quark 
counterpart (24). Heavy quarks lose most of their energy in the early stages of their 
trajectory when they are moving fastest. In contrast, from (45) we see that the light 
quark energy loss rate actually increases as x increases! This behavior, which was 
first pointed out in,®^ is reminiscent of a Bragg peak where the energy loss rate near 
thermalization is highest. Indeed, the energy loss rate (45) diverges as x —>■ Xstop-® 

Recently the light quark energy loss rate (45) was employed in a phenomenolog¬ 
ical hybrid model of jet quenching which included both strong and weak coupling 
physics. In the hybrid model the jet consists of a shower of partons which can 
split as the jet evolves. The splitting events were treated using perturbative QCD. 
Via soft interactions, in between splitings the partons can interact with the thermal 
medium and lose energy. The energy loss rate of the individual partons was treated 
using the light quark energy loss rate (45). The coefficient of the scaling in the 
stopping distance (26) was allowed to vary providing a one parameter fit to data. 
This gave a good fit to the data provided that the resulting stopping distances are 
a factor 3 to 4 times larger in QCD than that in SYM at the same temperature 
and energy. This is natural since SYM has more degrees of freedom for energetic 
partons to interact with than QCD. 

This will be expanded upon in a coming paper. 

® The rate that energy flows into hydrodynamic modes is only deflned over time and length scales 
^ 1/T. As such (45) should be regarded as the gradient expansion of the instantaneous jet energy 
loss rate. The gradient expansion of a function f(x) can diverge even if f(x) is itself finite. Hence, 
the instantaneous energy loss rate need not diverge as a; —> aj^top- 
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6. Discussion with future perspectives 

We hope to have shown some remarkable applications of holography to heavy ion 
collisions, while emphasizing that these results were obtained in a theory quite dif¬ 
ferent from QCD itself. Most notably almost all results are obtained in the limit 
of infinitely strong coupling. Nevertheless, many of the works presented promise to 
have a degree of universality: they may apply to a large range of strongly coupled 
gauge theories. One particular example of such a hope is figure 6, which could dis¬ 
play a universal rapidity profile of the plasma formed in high energy collisions in a 
strongly coupled gauge theories. Indications are that this profile can match exper¬ 
imental data for RHIC energies, while being perhaps too narrow at LHC energies. 

The latter may also be expected; it was already known that LHC collisions did 
not agree well with the holographic prediction of the total multiplicity of charged 
particles,as expanded upon in subsection 3.2. One natural explanation would 
be that at such high energies the coupling constant is intermediate, thereby inval¬ 
idating the assumption of an infinitely strong coupling. Indeed models with weak 
coupling tend to give wider rapidity distributions, so from RHIC to LHC we may 
be watching a cross-over. On the other hand it is necessary to make the AdS/CFT 
models more realistic before drawing a firm conclusion, which can for instance be 
done by having less symmetry, preferably containing event-by-event fluctuations in 
the initial profile. 

While this review naturally restricts to heavy ion collisions, we also wish to 
mention that holography is used for a much broader range of strongly coupled 
systems. These include non-fermi liquids^^'^ and holographic superconductors, 
whereby the latter are believed to possibly have some relation with high temperature 
superconductors. 

Lastly, we wish to express our excitement for the coming time, with several 
new experimental runs with different energies, different nuclei and higher statistics. 
Including features such as event-by-event distributions, non-trivial correlations, rare 
photon jets will give a very constraining data set, which will make it possible to 
distinguish all scenarios presented. This will hopefully give major lessons for both 
QCD, quark-gluon plasma physics, and applications of the gauge/gravity duality in 
a broad sense. 
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